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摘 要 : 以 无 核 白 葡萄 为 试 材 ,测定 了 其 在 不 同 季 节 (6 一 9 H). 不 同胞 间 CO, 浓 度 下 的 净 光 合 速 率 , 根据 贝 叶 
斯 方法 ， 结 合 蒙特 卡 罗 马 尔 科 夫 链 算法 对 光合 生化 模型 参数 进行 估算 ， 以 期 获得 不 同 季节 的 模型 参数 值 ， 并 
与 最 小 二 乘法 所 得 结果 对 比 , 探讨 贝 叶 斯 方法 在 解决 高 维度 复杂 模型 参数 估计 问题 中 的 可 行 性 和 葡萄 光合 作 


用 关键 参数 季节 变化 规律 。 2 
均 有 明显 的 季节 变化 特性 


吉 果 表明 ， 最 大 羧 化 速率 (Fnmaa)、 最 大 电子 传递 速率 (Jmaa)、 磷 酸 丙 糖 利用 速率 (TPU) 


， 出 现 先 增 后 减 的 趋势 ,8 月 达 最 高 ， 分 别 为 54.30 hmolm 2.S-1、88.45 umol-m ?s ! 


和 6.56 hmolm ?s 1!; 9 月 最 小 , 分别 为 34.66 umobkm ^s !.. 58.86 umol-m ^s ! a 4.38 humolm ^s ! ,. "FALE E 
(gm) 在 各 个 月 份 波动 不 大 , 6 一 9 月 分 别 为 5.16 humolm ?s Pa! 5.29 humolm ^s Pa! 5.39 umobm 7s Pa! 


和 5.41 numolm-?.s- Pa! . 


与 传统 的 最 小 二 乘法 相 比 ， 贝 叶 斯 方法 估算 的 Fnax 值 偏 小 , Jua TPU 和 gm 无 明显 


差异 。 同 时 贝 叶 斯 方法 估计 出 的 模型 参数 是 在 考虑 参数 先 验 信息 的 基础 上 获得 的 ， 生 化 意义 更 加 显著 。 试 验 
表明 ,光合 作用 生化 模型 (FVCB 模型 ) 在 应 用 于 光合 作用 模拟 时 ， 应 充分 考虑 其 参数 的 季节 变化 性 ; 结合 蒙特 
卡 罗 马 尔 科 夫 链 算法 的 贝 叶 斯 参数 估计 能 更 有 效 解决 FvCB 模型 中 参数 估计 问题 。 


关键 词 : 干旱 区 ; 葡萄 ; I 
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Abstract: The response of photosynthesis to CO» concentration can provide a number of important parameters related to envi- 


ronmental factors. Using white seedless grape as the tested material in this study, net photosynthetic rates of leaves were 


measured for different intercellular CO; concentrations during two typical growing seasons from June to September in 2014 
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and 2015. A widely used biochemical model (FvCB model) in the simulation of CO, and H5O gas exchange at the leaf scale 
was parameterized using data obtained from situ leaf-scale observations. In order to obtain the photosynthetic parameters val- 
ues, to explore seasonal variations in the photosynthetic parameters in different seasons and to discuss the feasibility and ad- 
vantage of the Bayesian method in solving high dimensional and complex model parameters estimation, the Bayesian approach 
was used to estimate the parameters of the FvCB model. In order to generate the Bayesian posterior probability distribution, a 
version of the Markov Chain Monte Carlo (MCMC) technique was used. In contrast, the least square procedure was used in the 
application of the same set of observational data. The results showed that maximum ribulose 1.5-bisphosphate carboxy- 
lase/oxygenase (Rubisco) carboxylation rate (V.max), potential light-saturated electron transport rate (Jmax) and the rate of use of 
triose-phosphates utilization (TPU) had evident seasonal variations which increased from June to August, and then decreased in Sep- 
tember. The maximum values were observed in August (54.30 umol-m ?s !, 88.45 umol-m ?^s ! and 6.56 umol:m ?s !, respectively) 


2.s and 4.38 umol m ?s ! 


and minimum values in September (34.66 umol:m ?^s !, 58.86 hmolm , respectively). The trend in mes- 
ophyll conductance (€m) was relatively stable in different months, with respective values of 5.16 hmolm ?s '-Pa'!, 5.29 
hmolm ^s '-Pa !, 5.39 umolm ^s Pa ', 5.41 umolm ^s '-Pa'! from June to September. In comparison with traditional least 
square method, the values of V4, estimated by the Bayesian method were relatively small and those of J4,,, TPU and gm had 
no obvious difference. Also because the estimated parameters by the Bayesian method were obtained after adequate considera- 
tion of prior information, each parameter was in biological sense obviously more meaning. As a consequence, it indicated that 
the Bayesian approach combined with Markov Chains and Monte Carlo (MCMC) sampling algorithm was an effective way of 
estimation of the parameters in the FvCB model. As the parameters in the FvCB model were different in different seasons, it 
was necessary to consider these variations in using the parameters in the FvCB model. 

Keywords: Arid region; Grape; Bayesian parameter estimation; Biochemical photosynthesis model; Photosynthetic parameter; 


Seasonal variation 


近年 来 ,光合 作用 的 模型 研究 得 到 人 们 日 益 关 
XU? HA, von Caemmerer 等 6-9 的 光合 作用 机 理 
模型 (以 下 简称 FvCB 模型 ) 因 具有 明确 生物 学 意义 
而 被 广泛 应 用 于 光合 作用 研究 中 ， 并 被 多 数 碳 循环 
模型 (如 SiB, CLM 等 ) 所 采用 。 在 该 模型 中 ， 最 大 羚 
化 速率 ( 矿 ss)、 最 大 电子 传递 速率 (ns)、 磷 酸 丙 糖 
利用 速率 (TPU) 是 表征 植物 光合 能 力 的 关键 参数 。 如 
何 利用 实测 净 光 合 速率 / 胞 间 CO, 浓度 变化 曲线 
(A/C; 曲线 ) 估 计 这 些 参数 不 仅 是 植物 生态 学 研究 的 
热点 ， 也 是 提高 陆地 碳 循环 模拟 精度 的 关键 [。 

前 人 以 不 同 植物 为 研究 对 象 对 模型 参数 估计 做 
TRZA. Harley tF] Wullschleger"! FH 4) Efi 
计 法 对 棉花 (Gossypium spp.) 和 109 种 C, 植物 的 
Visa, FEEDER E Ra 和 叶肉 导 度 gu 进行 了 估计 , 但 
此 类 方法 以 胞 间 CO; 浓度 分 压 等 于 20 Pa 为 节点 进 
行 数 据 的 经 验 性 分 割 ， 不 但 受 人 为 因素 扰动 较 大 ， 
而 且 当 数据 集 较 小 时 收敛 较为 困难 ; Dubois 等 (中 和 
Miao 等 上 结合 格 点 搜索 利用 非 线 性 最 小 二 乘法 同 
时 估计 出 Vans. Juss TPU 和 Ri， 但 该 方法 是 以 求 
误差 的 最 小 值 来 逼近 参数 真 值 ， 较 适用 于 连续 性 画 
数 的 参数 估计 ， 而 对 于 分 段 不 连续 且 高 维 的 FvCB 
模型 难以 获得 全 局 最 优 结果 ; Su 等 将 遗传 算法 应 
用 到 FvCB 模型 主要 参数 的 估计 当中 ,克服 了 传统 
友 代 方法 容易 陷入 局 部 最 优 解 的 缺陷 ,但 该 算法 后 
期 适应 度 容易 趋向 一 致 ， 使 优秀 的 个 体 在 产生 后 代 


时 优势 不 明显 ， 导 致 算法 进化 后 期 搜索 效率 较 低 ; 
近年 来 , Zhu 等 六 和 Feng 等 (1 先后 将 贝 叶 斯 方法 应 
用 到 光合 作用 模型 参数 估计 当中 ,， 较 其 他 方法 有 了 明 
显 的 优势 。 与 其 他 参数 估计 方法 相 比 ， 贝 叶 斯 方法 
在 充分 考虑 观测 误差 和 模型 结构 误差 的 基础 上 ， 获 
得 参数 的 分 布 范 围 ， 大 大 提高 了 拟 合 精度 。 然 而 将 
该 方法 应 用 于 光合 作用 生化 模型 参数 估计 ， 获 取 适 
用 于 我 国 西北 干旱 区 葡萄 的 模型 参数 的 研究 尚 鲜 见 
报道 。 

本 文 以 我 国 西北 干旱 绿洲 典型 农田 生态 系统 经 
济 作物 葡萄 (Vitis vinifera) 为 研究 对 象 ， 利用 便携 式 
光合 作用 -荧光 测量 系统 获取 不 同 季节 葡萄 净 光 合 
速率 随 胞 间 CO 浓度 变化 的 数据 。 在 此 基础 上 ,应 
用 贝 叶 斯 方法 对 FvCB 模型 参数 进行 估计 ， 并 分 析 
其 与 叶片 特性 的 关系 ,以 验证 贝 叶 斯 方法 在 解决 高 
维度 复杂 模型 参数 估计 问题 的 可 行 性 ， 并 估计 出 模 
型 参数 和 具有 的 生化 意义 ， 揭 示 葡 萄 光合 作用 关键 
参数 季节 变化 规律 及 其 与 叶片 特性 的 关系 ,增进 对 
干旱 绿洲 典型 经 济 作物 葡萄 光合 特征 的 认识 , 为 提 
高 农田 生态 系统 的 产量 提供 科学 指导 ， 推 动 贝 叶 斯 
方法 在 植被 光合 作用 模型 参数 估计 研究 上 的 应 用 。 


1 研究 区 概况 


试验 地 位 于 甘肃 省 敦煌 市 阳 关 镇 南湖 绿洲 西南 
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和 祁连山 、 北 塞 山 )， 西 接 塔克拉玛干 沙漠 ,观测 点 地 
理 坐 标 为 39°53'N, 94?07'E, 381 1 100~1 297 m; 深 
居 内 陆 ， 属 上 暧 温带 干旱 性 气候 , 日 照 时 间 长 , BA 


温差 大 ， 降 水 稀少 ， 蒸 发 量 大 ,年 均 降 水 量 36.9 mm, 


年 均 潜 在 蒸发 量 达 2486 mm", F ARAA 
3115-3 247 h， 年 均 气温 9.3 C; 土壤 主要 为 隐 域 性 
土壤 , 包括 沼泽 土 、 草 甸 土 和 盐 渍 土 等 。 观 测 点 位 
于 南湖 绿洲 西南 部 ， 属 农田 生态 系统 ,土壤 类 型 为 
绿洲 灌 耕 土 ， 和 葡萄 园 每 月 下 旬 进 行 一 次 人 工大 水 漫 
灌 , 确保 植物 在 充分 供水 条 件 下 生长 。 
2 材料 与 方法 
2.1 试验 材料 

试验 材料 为 ‘无 核 白 葡萄 *(Thompson Seedless), 
平均 冠 层 高 度 1.84m, 平均 胸径 3.2 cm, 树龄 13 a. 
随机 选取 4 棵 植株 , 利用 便携 式 光 合作 用 -荧光 测量 
系统 (GFS-3000) 进 行 CO; 响应 曲线 (4/C;) 测 定 。 观 测 
日 期 分 别 为 2014 年 的 6 月 、7 月 、8 月 、9 H, WA 
了 葡萄 的 开花 期 、 坐 果 期 、 浆 果 生 长 期 、 浆 果 成 熟 
期 4 个 典型 生长 期 。 观 测 时 间 处 于 每 月 中 下 旬 某 一 
青天 的 10:00 一 16:00。 观 测 叶 片 则 选取 位 于 冠 层 中 
上 部 的 成 熟 阳 生 叶片 ,用 于 不 同 季节 自然 条 件 下 的 
活体 测定 , 每 株 选 3 片 叶 子 进行 重复 试验 。 

测量 前 ,使 叶片 处 在 饱和 光照 条 件 下 
(PARiw=1 200 mol:m ^s 适应 30 min， 此 时 酶 活性 被 
完全 激活 。6、7、8、9 月 叶 室 温度 (7 分 别 设置 为 
25'C. 30 °C, 30 CH 20 C, 保证 了 与 当月 自然 环境 
温度 的 一 致 性 。 环 境 大 气压 强 (Pu 为 86 kPa, 叶片 温 
度 (Ti.n 设 定 范围 为 24~34 'C, 样品 室 相对 湿度 (RH) 控 
制 在 40%~65%, 均 处 在 当月 光合 作用 的 最 佳 范 围 内 ， 
气流 速率 为 750 hmols ,保证 了 CO 的 充分 吸收 。 在 
以 上 设置 不 变 的 情况 下 , 调整 CO, 绝 对 浓度 (COwaos) 来 
完成 不 同 植 株 、 不 同 叶 片 在 不 同 季节 的 A/C 曲线 测 
定 。COwavs 最 初 设 为 120 kPa, 而 后 为 100 kPa, 80 kPa, 
60 kPa, 50 kPa、 40 kPa、30 kPa、20 kPa、10 kPa 和 
5 kPa 共 10 个 水 平 ， 在 每 个 CO, 浓度 下 适应 2~3 min 
后 开始 测定 ， 且 每 个 梯度 记录 3 次 , 求 均值 。 
2.0 ”模型 简介 

FvCB 模型 将 整个 光合 作用 过 程 划分 为 3 个 限 
制 阶段 : 二 磷酸 核 酮 糖 羧 化 酶 限制 阶段 (Rubisco BR 
制 ) 二 磷酸 核 酮 糖 再 生 限 制 阶段 RuBP 限制 ) 和 磷酸 
丙 糖 利用 限制 阶段 (TPU 限制 )， 其 表达 式 为 FE-51: 


4n=min{4。， Ai, Ap? (1) 
AV ema (C/C +K.x(1+0/Ko)]-Ra — Q) 
A=Jx(Cf (4Cet8f )-Ra (3) 


^. 88H25 s 
Ay-3TPU-R4 (4) 
Ox. (ax Q--Jax)J ax QxJaax70 (5) 
C.-Ci-A/gm (6) 
R4-0.001 5 Vomax (7) 


式 中 : 4 为 净 光 合 速率 (umolm ^s ), A« A; 4, 43 
别 为 Rubisco 限制 、RuBP 限制 和 TPU 限制 阶段 的 
净 CO;[BSfbxk3E(umol:m ^s), Vi 为 最 大 羧 化 速 
X (umom ^s), J 为 电子 传递 速率 (humolm ^s), 
TPU 为 磷酸 两 糖 利 用 速率 (umolm-.s P), Ra 为 暗 呼 
吸 速 率 (umolm-s ,六 为 CO, 补偿 点 (Pa), 4 为 光合 
xkXE(umol-m ^s), KI. K。 分 别 为 羧 化 作用 和 加 和 氧 
作用 的 米 氏 常数 (kPa 和 Pa)，C. 和 o 分 别 为 Rubisco 
中 CO; 43 FE(Pa)fl O; 43 FE(21 kPa)! 9, Jna IRAE 
TE xÉ XE XE (umobm ^s), gmn 为 叶肉 导 度 
(umol:m ^s Pa ^), Ci 为 胞 间 CO» 7KI (Pa), 0 为 光 响 
应 曲线 斜率 (0.90), a 为 电子 转移 的 量子 产 率 (0.30), Q 
为 光合 作用 量子 通 量 密度 (0.093 hmolm 7s) 1, 
2.3 ”参数 估计 方法 

贝 叶 斯 定理 认为 ， 参数 的 后 验 分 布 与 先 验 分 布 
和 概率 密度 函数 的 乘积 成 正比 '"， 其 表达 式 为 : 


p(B1D) De pO | yx p(B) (8) 


式 中 : p 为 由 待 估 参 数 构成 的 参数 向 量 (包括 Vems 
Jus TPU 和 g,), DD 为 观测 数据 ，p(B81D) 为 后 验 分 
布 ，p(D|pB) 为 抽样 分 布 的 概率 密度 落 数 ，p(P) 为 
待 估 参数 6 的 先 验 概率 分 布 ，p(D) 为 随机 变量 的 
边缘 分 布 。 样 本 抽取 采用 Markov Chain Monte Carlo 
MCMO)Z; iE U?01, 

最 小 二 乘法 在 FvCB 模型 参数 估计 中 应 用 较 早 ， 
其 基本 思想 是 通过 最 小 化 误差 的 平方 和 寻找 数据 的 
最 佳 函 数 匹 配 00。 假 设 有 一 组 数据 x, y, i=1, 2，…， 
N, 且 已 知 该 组 数据 满足 了 (x) = fax), a 为 待定 参 
数 向 量 , 则 最 小 二 乘法 的 目标 是 求 出 一 组 待定 参数 
值 ， 可 以 定义 出 最 优化 问题 , 如 下 : 


N N 
J - min? Jy, -PE 2 min? -fla xn) (9) 
a ij 7o 


3 结果 与 分 析 
31 ”光合 作用 参数 估计 结 

以 无 核 白 葡萄 为 研究 对 象 ， 同 时 选用 最 小 二 乘 
法 和 贝 叶 斯 方法 对 FvCB 模型 运用 同一 套 观 测 数 据 
集 进行 参数 估计 (图 1)。 从 图 中 可 以 看 出 : PEEN 
叶 斯 方法 还 是 最 小 二 乘法 ， 参 数 季节 变化 显著 。 对 
于 参数 Vous Jeu 和 TPU， 其 值 呈 先 增 大 后 减 小 的 
趋势 ，8 月 达到 最 大 , 9 月 最 小 ; 而 gw" 则 是 6 月 和 9 
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月 略 大 于 7 月 和 8 月 ， 且 两 种 参数 估计 结果 的 季节 
变化 规律 一 致 。 值 得 注意 的 是 , 传统 的 最 小 二 乘法 获 
得 的 参数 估计 结果 是 一 个 定 值 ， 而 基于 贝 叶 斯 方法 的 
估计 结果 不 但 可 以 获得 参数 估计 结果 的 中 值 ,还 可 以 


"E 60 ? 
= fg 

£u 己 

长 g^ 

RO ŽŽ 


100 


Jains 
(umol:m™”-s™') 
oo 
© 


获得 参数 的 分 布 范围 。 此 外 ， 从 表 1 可 以 看 出 , 最 小 二 
乘法 估计 出 的 参数 Vi 在 4 个 月 中 均 大 于 贝 叶 斯 方法 
估计 结果 , 最 大 差别 可 达 15.51 umolm ^s (9 B), 而 
参数 JS. TPU 和 gm 无 此 现象 。 


e a 


* " 


磷酸 丙 糖 利用 速率 最 大 电子 传递 速率 


6 
人 
ze o * T 
S 4 
£ 
= 
Pe T T 
$6 
i 了 * $ 
Wb gs 4 
g'g 
= 5 e 最 小 二 乘法 
El Least square method 


6 7 


$ 


贝 叶 斯 方法 
Bayesian method 
8 9 


米 


月 份 Month 


1 不 同月 份 基于 贝 叶 斯 后 验 分 布 及 最 小 二 乘法 的 葡萄 光合 参数 估计 结果 (I 代表 贝 叶 斯 后 验 分 布 范围 ) 
Fig. 1 Posterior mean estimates results given by the Bayesian method and the least square method estimation results of grape pho- 
tosynthetic parameters in different months (I represents posterior distribution ranges based the Bayesian method) 
表 1 不 同月 份 基于 最 小 二 乘法 参数 估计 结果 和 基于 贝 叶 斯 方法 的 参数 先 验 信息 与 估计 结果 
Table 1 Parameters estimation results based on the least square method and the prior information and estimation results based on the 
Bayesian method in different months for grape photosynthesis 
最 小 二 乘法 贝 叶 斯 方法 Bayesian method 
"-— Least square method 先 验 信息 Prior information 后 验 分 布 Posterior distribution 后 验 中 值 Posterior median 
Parameter 月 份 Month TR ER 。 参考 文献 月 份 Month 
Lower | Upper 
6 7 8 9 limit limit Reference 6 7 8 9 6 7 8 9 
Vemex 36.65- 31.06- 32.46~ 23.76- 
(mom ^s 50.82 64.59 69.45 50.17 10 150 [7], [121 [13] 88.13 8847 88.98 90.69 48.21 52.98 54.30 34.66 
Jmax 55.57~ 57.41~ 63.89~ 47.26~ 
moim s 64.06 71.67 88.65 51.19 20 150 [7], [12],[13] 123.93 122.39 123.98 121.92 65.51 69.36 88.45 58.86 
TPU 3.30~  3.60-  428- 3.14~ 
(mol-m^«s^!) 4.84 5.06 6.40 3.70 0 9 [7], [12] 845 848 861 5.69 4.71 4.86 6.56 4.38 
£m 1.55». 1.465  1.53- 1.42~ 
(umol:m 2s ^p, !) 5.26 5.73 5.73 3.95 1 10 [7], [12] 9.40 9.36 9.32 9.34 5.601 5.29 5.39 5.41 
RMSE 0.58 1.39 1.12 0.75 — 1.06 0.70 1.26 0.76 


Vomax: 最 大 羧 化 速率 ; max: 


最 大 电子 传递 速率 ; TPU: 磷酸 丙 糖 利用 速率 ; gu HASE; RMSE: 3975 3801x296 Vom: maximum rate of 


Rubisco activity; Jmax: maximum electron transport rate; TPU: rate of triose-phosphates utilization; gm: mesophyll conductance; RMSE: root mean 


Square error. 


与 最 小 二 乘法 不 同 ， 贝 叶 斯 方法 是 一 种 参数 区 
闻 估 计 方 法 ， 其 目标 不 是 在 参数 可 行 区 间 内 找到 一 
个 使 模型 模拟 效果 最 好 的 最 优 参数 组 合 ， 而 是 在 充 
分 考虑 先 验 信息 误差 的 前 提 下 在 一 定 的 置信 水 平 上 


估计 模型 参数 的 分 布 区 间 。 从 不 同月 份 参数 的 先 验 
信息 及 后 验 分 布 结果 可 以 看 出 ( 表 1): 基于 先 验 信息 
的 贝 叶 斯 参数 估计 方法 可 以 有 效 地 缩小 给 定 的 参数 
先 验 分 布 范 围 。 定 义 UR(uncertainty reductions) 为 参 
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ZUR A XE EH A RSE (UR=1—Clposterior/CIprior Rr 
Clposterior 和 Clprior 分别 为 后 验 分 布 和 先 验 分 布 的 95% 
置信 区 间 )P0， 则 对 参数 Vau. Jma 和 TPU， 在 不 同 
月 份 其 不 确定 性 相对 减 小 量 近 50%， 最 高 达 7296; 
而 对 于 gm UR 值 相对 较 小 , 仅 略 大 于 10%( 表 2)。 由 
此 可 见 ， 参 数 Vous Jma 和 TPU KARIRA, 在 
贝 叶 斯 方法 下 可 被 很 好 地 估计 。 


R2 95% 和 置信 区 间 下 葡萄 光合 参数 不 同月 份 后 验 分 布 
相 较 于 先 验 分 布 的 不 确定 性 相对 减 小 量 

Table2 Relative uncertainty reductions in the length of 95% 

credible interval form prior to posterior distribution of grape 


photosynthesis in different months % 
e ud Vomax Jmax TPU gm 
6 63 47 43 13 
7 59 50 47 12 
8 60 54 52 13 
9 52 43 12 12 


V.max: 最 大 羧 化 速率 ; Ja 最 大 电子 传递 速率 ; TPU: 磷酸 丙 糖 
利用 速率 ; gm: HASE; RMSE: J97STRIRZÉo Voas: maximum rate 
of Rubisco activity; Jmax: maximum electron transport rate; TPU: rate of 


triose-phosphates utilization; gm: mesophyll conductance; RMSE: root 
mean square error. 


6 月 June 


3.2 ”基于 观测 值 和 模拟 值 的 模型 评价 

在 不 同胞 间 CO, 浓度 下 ， 基 于 贝 叶 斯 方法 的 净 
光合 速率 估计 值 CL4sm 和 净 光 合 速率 观测 值 dd,oo) 的 线 
性 回归 如 图 2 所 示 , 可 以 看 出 : 4 个 月 的 相关 系数 均 较 
高 , PHE 0.90 ME; 同时 , 各 月 观测 / 拟 合 趋势 线 趋 
近 于 1 : 1 直线 , 且 在 8 月 和 9 月 更 为 显著 ( 表 3)。 此 
外 , 与 最 小 二 乘法 相 比 ， 贝 叶 斯 方法 所 得 观测 / 拟 合 数 
据点 整体 分 布 在 更 加 靠近 1 : 1 直线 的 区 域 ， 充分 表明 
贝 叶 斯 方法 县 有 更 强 的 收敛 性 ， 获 得 了 更 高 精度 的 参 
数 估计 结果 。 

光合 作用 对 胞 间 CO, 浓 度 的 响应 曲线 (4/C; 曲线 ) 
是 分 析 光 合作 用 机 理 的 重要 指标 。 将 不 同方 法 的 参 
数 估 计 结 果 应 用 于 FvCB 模型 , 4/C; 曲线 拟 合 结果 如 
3 所 示 。 从 图 中 可 以 看 出 : 对 于 贝 叶 斯 方法 的 拟 
合 结 果 ， 净 光合 速率 (4,) 随 胞 间 CO; 浓度 (Cj) 的 升 高 
呈现 出 明显 的 3 个 阶段 , 与 FvCB 模型 结构 相 了 吻合 ， 
B]: 当 C; 相 对 较 低 时 ，4, 随 Ci 的 升 高 呈 直 线 快 速 上 
升 ， 此 时 处 于 光合 作用 第 1 阶段 一 一 Rubisco 限制 阶 
段 ; 随 着 Ci 的 变 大 , 如 随 Ci 升 高 的 增 速 逐渐 变 小 ,此 
时 处 于 第 2 阶段 一 一 RuBP 限制 阶段 ; 当 Ci; 继续 增 大 


i or 
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[s 

$ o 
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净 光 合 速率 拟 合 值 4sn 
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0 L L 
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净 光 合 速 率 观测 值 AL, (umol m ^s") 


贝 叶 斯 方法 


Bayesian method 


o 最 小 二 乘法 1: 


Least square method 


1 直线 


1: 1 line 


2 不 同月 份 葡萄 净 光 合 速 率 观测 值 和 贝 叶 斯 方法 及 最 小 二 乘法 拟 合 的 线性 回归 
Fig.2 Regressions of grape net photosynthesis rate between measured value (Anovs) and modeled value (Ansim) by the Bayesian 
method and the least square method in different months 
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表 3 基于 贝 叶 斯 方法 的 葡萄 不 同月 份 净 光合 速率 观测 值 (4wovs) 和 拟 合 值 (4wsim) 的 统计 信息 
Table3 Statistical information of measured (Anops) and estimated (Ansim) values of net photosynthesis rate of grape based on the 
Bayesian method in different months 
m 2 
月 从 回归 方程 平均 偏差 纳什 系数 (NSE) — 归 一 化 平均 绝对 误差 (NMAHB) 
Month Regression equation ub ci 2 Mean bias Nash coefficient Normalized mean absolute error 
6 dnsim=0.834nobs+0.79 0.99 1.04 0.95 0.12 
7 dnsim=0.894nobs+0.90 0.99 0.65 0.98 0.07 
8 Ansim=0.92Anobs+0.86 0.99 1.00 0.96 0.08 
9 Ansim=0.92Anobs+1.26 0.99 0.59 0.97 0.07 
6 月 June 7 月 July 8 月 August 9 月 September 
25r 25r 
2 20} zx 
A 
- 
I5F 15 上 P 
p pea ei 
个 10r 10F : 
"LA 4 l E T V ， 
5 | 多” 贝 叶 斯 方法 "E 贝 叶 斯 方法 贝 叶 斯 方法 uum 贝 叶 斯 方法 
El or* Bayesian method Bayesian method 0r Bayesian method Bayesian method 
Ş 25 30 25 
x : / 7 
ES 20r I / 20|- / 
$ iss-_/ A oct ee ete ds] 7 
ES -- 
^m qr 10|- 
. e 10 Ins . 
5p 最 小 二 乘法 5| /最 小 二 乘法 ° Kuh Ri 5- V^ 最 小 二 乘法 
L Least square Least square east square - Least square 
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Observed value 


Rubisco-limited stage 


RuPB-limited stage TPU-limited stage 


3 不 同月 份 葡萄 净 光 合 速率 随 细胞 间 二 氧化 碳 浓度 (CD 变化 趋势 图 


Fig.3 Prediction of net photosynthesis rate (4n) of grape as a function of intercellular carbon dioxide concentration (Ci) in different months 


到 一 定数 值 时 ，4, 不 再 随 C; 的 升 高 而 变 大 ， 而 是 趋 
于 稳定 ， 此 时 处 于 第 3 阶段 一 一 TPU 限制 阶段 。 且 4 
个 月 第 1 阶段 和 第 2 阶段 转折 点 对 应 的 Ci 分 别 为 
31.68 Pa、48.19 Pa, 50.82 Pa 和 35.67 Pa， 出 现 了 明 
显 的 季节 差异 性 。 此 外 , 不 同月 份 最 大 净 光 合 速率 
分 别 为 : 13.49 umol-m ?:s 15.32 umol:m ?^:s 4, 18.49 
umol:m ^:s ^, 12.28 umo m ^s 1， 与 参数 估计 结果 有 
相同 的 季节 变化 规律 。 而 对 于 最 小 二 乘法 的 拟 合 结 
果 ， 尽 管 确 保 了 均 方 根 误差 保持 在 较 小 的 范围 ， 由 
于 估计 时 参数 之 间 的 相互 作用 ， 只 拟 合 出 了 二 磷酸 
核 酮 糖 再 生 限制 阶段 (RuBP 限制 ) 和 磷酸 丙 糖 利用 限 
制 阶段 (TPU 限制 )， 导 致 了 生理 生态 学 意义 不 够 明 
显 。 可 见 ， 相 较 于 最 小 二 乘法 ， 贝 叶 斯 方法 可 以 估计 
出 具有 生化 意义 的 模型 参数 ,参数 优化 后 的 FvCB 模 
型 可 用 于 模拟 我 国 西北 干旱 绿洲 作物 的 光合 速率 。 


4 讨论 
无 核 白 葡萄 为 多 年 生 木 本 ,参数 Vosaxss Jmax 和 


TPU 的 估计 范围 与 Wullschleger 结果 相似 ，Wnull- 
schleger 选取 了 109 种 C, 植物 ， 基 于 A/C 曲线 的 非 
线性 回归 技术 估计 了 不 同 植被 类 型 的 参数 值 , 给 出 
参数 变化 范围 : Vs 为 6~94 hmolm ^s 1!，Jiwax 为 
17-372 pmol:m“s !, TPU 为 4.9-20.1 jmol-m™s !, 
其 中 一 年 生 草 本 植物 参数 值 大 于 多 年 生 木 本 植物 。 

本 文 利 用 贝 叶 斯 方法 精准 地 估计 了 模型 参数 ， 且 结 
果 在 前 人 所 给 结果 范围 内 。 同 时 ,传统 参数 估计 方 
ih: 如 Harley 等 名 利用 分 段 方法 估计 出 了 Vema Ra 
和 gm， 但 该 方法 将 A/C; 曲线 中 Rubisico 限制 阶段 和 
RuBP 限制 阶段 的 转折 点 设 为 定 值 , 而 本 文 贝 叶 斯 
拟 合 结果 显示 出 转折 点 明显 的 季节 差异 性 ， 反映 出 
分 段 估计 法 受 人 为 分 割 的 扰动 会 产生 较 大 误差 ; 
Sharkey 等 {用 微软 电子 表格 同时 估计 出 Vsus. Jas 
TPU, Rafi gm; Dubois 等 00 和 Miao 等 5 结合 了 格 
点 搜索 和 非 线 性 最 小 二 乘法 估计 出 Vases Juss TPU 
和 R4， 这 些 方法 都 是 派生 的 最 优化 方法 , 不 但 计算 
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量 较 大 、 对 初始 值 的 设 定 十 分 敏感 ， 而 且 适 用 于 连 
续 画 数 的 参数 估计 ,而 对 于 分 段 不 连续 的 FvCB 模 
型 难以 得 出 全 局 最 优 结果 。 相 反 ， 贝 叶 斯 方法 不 需 
要 通过 求 导 函数 实现 目标 函数 的 最 小 化 ， 同 时 又 结 
合 了 先 验 信 息 , 所 以 该 方法 在 解决 高 维 复杂 不 连续 
函数 的 参数 估计 间 题 上 有 很 大 优越 性 。 本 文 首次 将 
该 方法 应 用 在 我 国 西部 极端 干旱 的 敦煌 绿洲 ， 准 确 
估计 了 当地 主要 经 济 作物 葡萄 在 典型 生长 季 (6 月 、7 
H. 8 月 和 9 月 ) 内 不 同 阶段 (开花 期 、 坐 果 期 、 浆 果 
生长 期 、 浆 果 成 熟 期 ) 的 模型 参数 ， 对 提高 农田 生态 
系统 的 产量 有 一 定 指导 , 为 陆地 生态 系统 碳 循环 模 
型 研究 提供 一 定 的 借鉴 。 

不 同月 份 参数 估计 结果 表现 出 明显 的 季节 差异 ， 
与 葡萄 在 不 同 季节 氮 分 配 策略 相关 。 无 核 白 葡萄 为 
C; 植 物 , 约 60%~80% 的 叶片 氮 以 核酸 和 酶 的 形式 存 
在 于 叶片 中 ?小 。 前 人 研究 表明 ,叶片 含 氮 量 与 光 
合作 用 能 力 存 在 正 相 关 关 系 “1。 在 光合 能 力 较 强 
的 月 份 , 所 需 氮 素 多 ， 酶 活性 强 ; 反之 , 光合 能 力 较 
弱 的 月 份 ,， 氮 素 消耗 较 少 ， 酶 活性 较 弱 ?1。 Vou. 
是 Rubisco 限制 阶段 的 最 重要 参数 ， 此 时 CO, 的 固定 
受 Rubisco 活性 的 制约 ,所 以 在 8 月 ,参数 V... 也 
处 于 最 大 值 ; 而 在 以 9 月 份 为 代表 的 浆果 成 熟 期 ， 
叶片 趋 于 衰老 ， 酶 活性 较 弱 ， 在 无 大 量 气 素 消 耗 时 ， 
叶片 氮 向 木质 部 转移 贮存 ， 参 数 Vau, 也 出 现 较 低 
值 。 同 理 , 分 别 作 用 于 RuBP 限制 阶段 和 TPU 限制 
阶段 的 参数 Jma 和 TPU 也 表现 出 相同 的 季节 变化 特 
性 ,而 参数 gn 则 受 温度 和 太阳 辐射 影响 较 大 , 7 月 和 
8 月 相 较 于 6 月 和 9 月 具有 较 高 的 温度 和 较 强 的 太 
阳 辐 射 ， 此 时 叶片 为 防止 大 量 蒸发 失 水 ,会 出 现 气 
孔 关 闭 情况 ， 从 而 导致 较 小 的 气孔 导 度 。 


5 结论 


1) 基 于 贝 叶 斯 方法 估计 出 FvCB 模型 中 4 个 主要 参 
数 范围 : 最 大 其 化 速率 (Vs) 为 23.76~90.69 pmolm*s 
最 大 电子 传递 速率 (i,) 为 47.26~123.98 umoblm ^s !, W 
酸 丙 糖 利用 速率 (TPU) 为 3.14~8.61 umol:m ^s !, 叶肉 导 
度 (@) 为 1.42-9.40 umol-m ?s -Pa 

2) 参 数 大 小 存在 明显 的 季节 变化 规律 ， 最 大 羚 
REV) RA BFE ER E (Jna) M R A 
利用 速率 (TPU) 参 数 大 小 出 现 先 增 后 减 的 趋势 ,8 月 
最 大 9 月 最 小 , 气孔 导 度 (gm 则 6 月 和 9 月 略 大 于 7 
月 和 8 月 。 

3) 与 传统 的 最 小 二 乘法 相 比 ， 贝 叶 斯 参数 估计 
能 够 有 效 解决 高 维度 复杂 不 连续 模型 的 参数 估计 问 
题 ， 可 以 基于 观测 数据 有 效 估计 出 上 县 有 生化 意义 的 


合理 参数 值 。 
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